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Résumé — On aborde la régulation du niveau d'eau dans un 
aménagement hydraulique, soumis à de fortes contraintes, 
par la commande sans modèle. Les nombreuses simulations nu- 
mériques fournissent d'excellents résultats, obtenus grâce à 
des algorithmes robustes et simples. 

Abstract — A new concept called Model-Free Control is applied 
to hydroelectric run-of-the river power plants, with severe 
constraints and operating conditions. Numerous computer 
simulations display excellent results, which are obtained 
thanks to simple and robust algorithms. 

Mots-clés — Biefs, canaux, surface libre, commande sans mo- 
dèle, PID intelligents. 

Keywords — Reaches, channels, free surface, model-free 
control, intelligent PID controllers. 



I. Introduction 

La régulation du niveau d'eau dans un aménagement hy- 
draulique (barrages, retenues, rivières, canaux, . . .), impor- 
tante pour la production d'électricité, la navigation, l'irri- 
gation, et bien d'autres usages, a reçu une attention consi- 
dérable, selon des points de vue les plus variés, dans la 
littérature automaticienne récente : voir, par exemple, [3J, 

m, mi, m, m, 0, uni, im, en, uni, eu, eu, eé, hh, 

[30], et les récents tours d'horizon [TÏÏJ, [25], et [33J. La 
difficulté tient à la nature même du procédé : 

- c'est un système hydraulique d'écoulement à surface 
libre, enchaîné (voir Fig. Q]) ou non, avec géométrie 
quelconque du volume d'eau ; 

- sa taille est grande (plusieurs kilomètres de long) ; 

- il est non linéaire, multivariable et distribué dans 
l'espace (équation aux dérivées partielles de Saint- 
Venant) ; 

- il existes des perturbations aléatoires dont l'amplitude 
peut être importante, mais inconnue, dues à des ap- 
ports ou prélèvements d'eau. 

L'embarras croît encore si 

- le niveau d'eau à réguler est éloigné de l'organe ré- 
glant (une vanne, par exemple), compte-tenu des re- 
tards variables mais aussi de phénomènes hydrauliques 
complexes de basculement des plans d'eau, en fonction 
du débit ; 

- le cahier des charges exige une robustesse des perfor- 



mances sur toute la plage de fonctionnement et vis- 
à-vis des variations des caractéristiques de l'aménage- 
ment (temps de retard variables, superficie du réser- 
voir variable en fonction du niveau et de l'envasement, 
usure des organes réglants, . . .). 
Cette communication, qui complète une précédente |23) . 
traitant d'un cas simplifié, propose pour des aménagements 
hydroélectriques en cascade, séparés par des biefs (voir Fig. 
Q]), la nouvelle « commande sans modèle » ([T5], [13]), facile 
à mettre en œuvre, et éminemment robuste, dont nous rap- 
pelons l'essentiel au § [TTJ Après une description du disposi- 
tif et de sa régulation au S Mil le S IIVI fournit d'excellentes 
simulations numériques, respectant plusieurs scénarios. Le 
§ [V] conclut au grand intérêt pratique de notre approche, 
tout en posant la question du rôle de la dimension infinie. 

Remerciements. Travail sous l'égide de deux contrats, intitulés 
« Commande sans modèle pour aménagements hydroélectriques en 
cascade », entre, d'une part, EDF et, d'autre part, l'Ecole polytech- 
nique, 1TNRIA et le CNRS. Le système de régulation, qui y est décrit, 
a fait l'objet d'une demande de brevet (n° FR0858532), déposée par 
EDF et l'École polytechnique le 12 décembre 2008. 

IL Rappels 

A. Commande sans modèle et i-PID 

La commande sans modèle repose sur une modélisation 
locale, sans cesse réactualisée, à partir de la seule connais- 
sance du comportement entrée-sortie. A l'équation diffé- 
rentielle inconnue, linéaire ou non, 
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(i) 



décrivant approximativement le comportement entrée- 
sortie, on substitue le modèle « phénoménologique », va- 
lable sur un court laps de temps, 



y 



F + /3u 



ou 



- l'ordre de dérivation v, en général 1 ou 2, donc différent 
de l'ordre de dérivation b de y en )Q, est fixé par 
l'opérateur ; 




Fig. 1: Schéma d'un bief 



- le paramètre constant (3 est fixé par l'opérateur de 
sorte qu'en ([T]) fiu et y( v > aient même ordre de gran- 
deur. 

La valeur de F à chaque instant se déduit de celles de u et 
de y^ v \ obtenus par dérivateurs numériques. On obtient le 
comportement désiré, si, par exemple, v = 2 en ([2]), grâce 
au correcteur PID intelligent, ou, en abrégé, i-PID, 



K D è 



(3) 



ou 



- y* est la trajectoire de référence de la sortie, obtenue 
selon les préceptes de la commande par platitude ; 

- e = y — y* est l'erreur de poursuite ; 

- Kp, Ki, Ko sont les gains de réglage. 

Remarque 1 II est aisé de déterminer ces gains car, avec 
©, © se ramène, contrairement aux PID classiques, à un 
intégrateur pur du second ordre. 

Remarque 2 Voir 

- J12]/ . \1S$ pour plus de détails, 

- pour les liens avec les PID usuels, 

- U!, M, fZSF, JZ2P, HE E22F Pour plusieurs autres illus- 
trations concrètes. 

B. Dérivateurs numériques 

B.l Un calcul simple 

Avec les notations classiques du calcul opérationnel {cf. 
[32]), il correspond à p(t) = a<j + ait, ao,a\ G R, pour 
t > : 

Pis) = ^ + ^ (4) 
s s 

On cherche à éliminer car on veut estimer ai. Pour cela 
on multiplie Q par s : 



sP(s) = a + 



ai 



puis on dérive l'expression obtenue par rapport à s pour 
supprimer ag : 

dP(s) _ ai 
ds s 2 



P(s) 



Avant de revenir au domaine temporel, une multiplication 
par s~ N , avec N > 1, N = 2 par exemple, est nécessaire 
pour obtenir uniquement des intégrales : 



>P(s) + s-^ = -s- 
as 



ai. 



On revient au domaine temporel en rappelant (cf. |32| ) que 
correspond à la multiplication par —t : 

6 / Tx(r)àT — j j x(n)AHiàT 

\Jt JtaJta J 



ai = 



6 ( / TX ^ dT ~ / (* ~ T M T ) dr ^ (5) 
? 

6 / (tx(t) — (t — r)cc(r))dr 

Jtg 

t 3 



B.2 Filtres dérivateurs 

Les propriétés du calcul opérationnel permettent de gé- 
néraliser §5§ en écrivant tout estimateur algébrique d'une 
dérivée d'ordre quelconque d'un signal x(t), analytique au- 
tour de 0, sous la forme d'un filtre J t " vo(t)x(t)àt, où vo(t) 
est un polynôme temporel caractérisant ledit estimateur. 
Les calculs du § III-B.1I permettent de comprendre la va- 
riété possible d'estimateurs que l'on peut obtenir par cette 
méthode, en modifiant, par exemple, l'ordre de troncature 
du développement de Taylor. 



Remarque 3 Pour compléter ce rapide résumé, inspiré de 
\15y . renvoyons à \2^j pour des développements théoriques 
importants, y compris pour l'implantation algorithmique. 
Voir \1J$ pour plus de références et, surtout, pour les pro- 
grès substantiels que permet cette dérivation en automa- 
tique non linéaire, avec modèle connu, aux incertitudes pa- 
ramétriques près. 




C. Commande sans modèle 
Nous mettons ici en œuvre ^ , avec v — 1 , sous la forme 



Fig. 2: Schéma bloc du modèle 

III. Description du dispositif 

A. Généralités 

Comme l'indiquent les schémas [T] et [2j on cherche à maî- 
triser le niveau d'eau z r du bief i, entre les usines i — 1 et 
i, avec 

- commande u, ou débit de sortie Q s , échantillonnée- 
bloquée et saturée en position et vitesse ; 

- deux entrées exogènes Q e et W, perturbatrices ; 

- insensibilité du capteur de niveau qui a une résolution 
de 1 cm. 

La Fig. fournit le schéma de commande : 
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Bief 
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Fig. 3: Schéma bloc de la commande désirée 



B. Traitement du retard 

Le débit Q e modifie le niveau z r . loin de l'usine i, avec 
un retard mal connu et, qui plus est, variable, donc diffi- 
cile, sinon impossible, à maîtriser avec les outils théoriques 
actuels d'identification et de commande. On y pallie en re- 
construisant un niveau z a à partir de mesures proches de 
l'usine i — I , par deux lois possibles : 

f . l'une repose sur la connaissance de niveaux 



f{Ze) 



2. l'autre sur celle du débit entrant 



g(Q e 



fournies par EDF et qui résultent de lois empiriques. 

La trajectoire de référence, intermédiaire, z*, détermi- 
née selon les principes de la platitude, est complétée par 
un correcteur piloté par l'erreur de poursuite sur z r . Nous 
réalisons, ainsi, une replanification en ligne automatique de 
la trajectoire de référence. Cette trajectoire corrigée z* de- 
vient la consigne d'une boucle interne, à commande sans 
modèle. 



ay = F — u 



(0) 



qui a l'avantage de donner un sens clair de débit entrant 
à F, information importante pour l'exploitation de la cen- 
trale hydraulique, en faisant apparaître la différence F — u, 
où u est le débit sortant. La boucle est fermée, d'après ©, 
par 

u = F - ai/* + oLT(e) 

où PI(e) = Kpe + Kj J e est un correcteur PI, piloté par 

e = y - y*- 

IV. Quelques résultats de simulation 

A. Scénarios 

Trois scénarios, proposés par EDF sur l'aménagement 
hydroélectrique de Fessenheim, permettent d'évaluer les 
performances de notre commande. Cet aménagement au 
fil de l'eau (capable de turbiner 175, 5 MW en pointe) fait 
partie des 10 centrales EDF (avec sa filiale EnBW) enchaî- 
nées en cascade sur le Rhin et régulées en niveau. Il est 
caractérisé par un bief de 15 kilomètres de long et d'une 
hauteur de chute d'environ 15 mètres. Le principal objectif 
à satisfaire est le respect de la contrainte 

z* - 10 cm < z r < z* + 10 cm (7) 

en un point distant de 7,5 kilomètres de l'actionneur. 

1 . Le premier scénario représente la fin d'une crue avec de 
fortes variations de débit. 

2. Le second, plus doux, correspond à une situation nor- 
male plus fréquente ; 

3. Quant au troisième, il est académique afin de mettre en 
lumière les comportements lors de saturations. 

Les scénarios 1 et 2, qui sont tout à fait réalistes, ont une 
durée de 4 jours. Les figures peuvent donner l'illusion, donc, 
d'une commande très dynamique, ce qui n'est pas vraiment 
le cas. 

Les perturbations sont constituées, d'une part, par un 
biais égal à 0.03Q e + 10 pour prendre en compte les erreurs 
de débit (consigne - mesure) et, d'autre part, par des sassées 
de 100 m 3 /s, reproduisant le comportement d'une écluse. 
Ces dernières sont très violentes puisqu'elles ne durent que 
15 minutes, c'est-à-dire 7 échantillons. 

B. Réglage de la boucle externe 

Nous avons rencontré des difficultés à régler le correcteur 
de la boucle externe car seule une méthode empirique est 
applicable. Nous avons respecté la règle classique stipulant 
que la dynamique du correcteur de la boucle externe est 
moins rapide que la boucle interne. 

C. Commande échantillonnée bloquée à 2 minutes 

Les Figures [4] et [5] présentent les résultats obtenus dans 
le cas où la commande est bloquée toutes les 2 minutes, en 
respectant les différents scénarios. Aucune contrainte n'est 
violée : le contrat est rempli ! 



Remarque 4 Nous avons repris ces mêmes simulations en 
augmentant la fréquence d 'échantillonnage, en passant à 1 
minute. On améliore considérablement ([7]) en obtenant : 
z* — 5 cm < z r < z* + 5 cm. 

D. Anti- emballement 

Un anti-emballement, ou anti-windup, est mis en place 
sur chaque boucle de régulation (boucles interne et ex- 
terne). C'est pourquoi, sur la Figure [51 la commande et 
la replanification de z* décroissent immédiatement après 
la diminution du débit d'entrée Q e . 

E. Consigne anticipée 

Il est, bien sûr, possible d'anticiper la consigne de z r : 
ceci revient à prévoir les changements de consigne et à les 
appliquer, afin de considérer le temps séparant l'évolution 
de z de son effet sur z r . Afin de ne ne pas masquer les résul- 
tats de la régulation proposée, cette idée n'est pas mise en 
œuvre ici. De fait, on voit apparaître, dans la commande à 
2 minutes, un retard d'environ 45 minutes que nous pour- 
rions largement atténuer. 

V. Conclusion 
A . Intérêt pratique de la commande sans modèle 

La stratégie proposée possède de nombreux avantages : 

- Le retard, cause d'instabilité dans bien des lois de com- 
mande, est écarté grâce à des redondances analytiques 
approximatives, déduites de la physique du bief. Les 
erreurs, qui en résultent, sont corrigées par une boucle 
externe de régulation. Une nouvelle trajectoire de ré- 
férence est calculée à chaque instant pour le niveau 
aval. 

- On tire toujours profit de la grande réactivité de la 
commande sans modèle dans la boucle interne, qui as- 
sure ainsi la poursuite de la trajectoire sans connaître 
le modèle du bief. 

Par ses propriétés de robustesse, d'adaptabilité et de sim- 
plicité, la commande sans modèle apporte des performances 
remarquables avec un temps de mise au point très court, 
comparé aux systèmes avancés de commande, que ce soit 
en étude de simulation ou sur site. Elle semble donc parti- 
culièrement adaptée au milieu industriel. 

L'algorithme de commande proposé est innovant dans le 
domaine de la régulation de niveau des canaux découverts. 
Les performances obtenues en poursuite, en rejet de pertur- 
bations, et en robustesse sont remarquables, compte tenu 
de la sévérité des scénarios de simulation appliqués. Le fait 
de maintenir le niveau au point milieu du bief - soit à 7,5 
kilomètres de l'actionneur - dans une bande de ±10 cen- 
timètres, avec des perturbations inconnues, relève en effet 
d'une vraie prouesse. La commande sans modèle appliquée 
à la régulation de niveau constitue une régulation indus- 
trielle non seulement en raison des performances atteintes 
mais aussi, et surtout, pour ses qualités intrinsèques : 

- facilité de mise en oeuvre (structure de type PID feed- 
forward) , 

- faible sollicitation de la charge du calculateur, 

- temps de mise au point réduit, grâce, notamment, au 
faible nombre de paramètres à régler[U, 

1. 5 si l'on fait la somme des paramètres dans J2J et J3}. Si l'on part 



- auto-adaptation par rapport aux variations du pro- 
cédé, 

- maintenance aisée, grâce à un algorithme très simple. 
L'intérêt pour le producteur hydraulique EDF est de dispo- 
ser d'une commande capable de réguler un niveau proche 
ou distant du vannage (organe réglant), tout en garantis- 
sant les performances désirées sur toute la plage de fonc- 
tionnement de l'usine hydroélectrique, malgré la présence 
de perturbations imprévisibles. Ces atouts permettront cer- 
tainement de réduire les coûts des projets de construction 
ou de rénovation d'installations hydrauliques, en diminuant 
la durée des essais de mise en service. 

B. Quel rôle pour la dimension infinie ? 

Des techniques élémentaires suffisent, ici, pour traiter 
d'un sujet qui relève de la dimension infinie (systèmes 
à retards et régis par des équations aux dérivées par- 
tielles). Reste à savoir si cette constatation peut s'étendre 
à d'autres applications dont la modélisation naturelle re- 
quiert a priori la dimension infinie. C'est, sans aucun doute, 
une question méthodologique, voire épistémologique, fon- 
damentale pour le futur de l'automatique et, peut-être, de 
plusieurs autres domaines d'ingénierie et de mathématiques 
appliquées. Seul l'examen de multiples cas concrets, posés 
par les praticiens, permettra d'y répondre. 
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Fig. 4: Scénarion 1 avec poursuite et rejet de fortes perturbations : période de 2min 
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Fig. 5: Scénarion 2 avec poursuite et rejet de perturbations assez lentes : période de 2min 
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Fig. 6: Scénarion 1 avec poursuite et rejet de perturbations avec actionneur saturé : période de 2min 
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